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The Photonic Force Microscope (PFM) is an opto-mechanical technique based on an optical trap 
that can be assumed to probe forces in microscopic systems. This technique has been used to measure 
forces in the range of pico- and femto- Newton, assessing the mechanical properties of biomolecules 
as well as of other microscopic systems. For a correct use of the PFM, the force field to measure 
has to be invariable (homogeneous) on the scale of the Brownian motion of the trapped probe. 
This condition implicates that the force field must be conservative, excluding the possibility of a 
rotational component. However, there are cases where these assumptions are not fulfilled Here, we 
show how to improve the PFM technique in order to be able to deal with these cases. We introduce 
the theory of this enhanced PFM and we propose a concrete analysis workflow to reconstruct the 
force field from the experimental time-series of the probe position. Furthermore, we experimentally 
verify some particularly important cases, namely the case of a conservative or rotational force-field. 

PACS numbers: 05.40.Jc, 87.80. Cc, 07.10. Pz, 47.61.-k 

Keywords: Brownian motion, Photonic Force Microscope, Singular points, Stability analysis 

I. INTRODUCTION 

A focused optical beam - an optical tweezers - permits one to manipulate a wide range of particles - including atoms, 
molecules, DNA fragments, living biological cells, and organelles within them - with applications to many areas - such 
as molecular biophysics, genetic manipulation, micro-assembly, and micro-machines [US IS]- One of the most exciting 
applications has been the possibility to investigate and engineer the mechanical properties of microscopic systems - 
using, for example, optical traps as force transducers for mechanical measurements in biological systems [H [5j |6j [7j [8] . 

In the early 90s various kinds of scanning probe microscopy were already established. The Scanning Tunneling 
Microscope (STM) |9J permits one to resolve at the atomic level crystallographic structures [TU] and organic molecules 
[TT] . The Atomic Force Microscope (AFM) [T^] has been successfully employed to study biological and nano-fabricated 
structures, overcoming the diffraction limit of optical microscopes. Furthermore, they developed from pure imaging 
tools into more general manipulation and measuring tools on the level of single atoms or molecules. However, all 
these techniques required a macroscopic mechanical device to guide the probe. 

A new kind of scanning force microscope using an optically trapped dielectric microsphere as a probe was proposed 
in |T3l [14]. This technique was later called Photonic Force Microscope (PFM) [15]. In a typical setup, the probe is 
held in an optical trap, where it performs random movements due to its thermal energy. The analysis of this thermal 
motion provides information about the local forces acting on the probe. The three-dimensional probe position can 
be recorded through different techniques which detect its forward or backward scattered light. The most commonly 
used are a quadrant photodiode, a position sensing detector, or a high-speed video-camera [3]. The PFM provides the 
capability of measuring forces in the range from femto- to pico-Newton. These values are well below those achieved 
with techniques based on micro-fabricated mechanical cantilevers, such as AFM |16j . 

For small displacements of the probe from the center of an optical trap, the restoring force is proportional to the 
displacement. Hence, an optical trap acts on the probe like a Hookeian spring with a fixed stiffness, which can be 
characterized with various methods [3J [T7] . The correlation or power spectrum method, in particular, is considered 
the most reliable [17 , allowing one to determine the trap stiffness by applying Boltzmann statistics to the position 
fluctuations of the probe, relying only on the knowledge of the temperature and the viscosity of the surrounding 
medium m m EH HI US] ■ 

Assuming a very low Reynolds number regime HI] , the Brownian motion of the probe in the optical trap is 
described by a set of Langevin equations: 

7 r(i) + Kr(f) = %/2\D 7 h(t) , (1) 

where r(t) = [x(t),y(t), z(t)} is the probe position, 7 = GttRj] is its friction coefficient, R is the probe radius, 77 is the 
medium viscosity, K is the restoring force matrix, ^/2Djh(t) — ^/2D r )[h x (t),h y (t),h z (t)} is a vector of independent 
white Gaussian random processes describing the Brownian forces, D = keT/j is the diffusion coefficient, T is the 
absolute temperature, and fcs is the Boltzmann constant. The orientation of the coordinate system can be chosen in 
such a way that the restoring force is independent in the three directions, i.e. K = diag(fe x , k y , k x ). In this reference 
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FIG. 1: (Color online) Examples of physical systems that produce forces fields that can not be correctly probed with a classical 
PFM because they vary on the scale of the Brownian motion of the trapped probe, i.e. in a range of 10s to 100s of nanometers 
: (a) forces produced by a surface plasmon polariton in the presence of a patterned surface on a 50 nm radius dielectric particle 
(reproduced from [3T]); (b) electromagnetic near-field of a 10 nm wide gold tip in water illuminated by a 810 nm monochromatic 
light beam (reproduced from 32 ); and (c) vector force field acting on a 500 nm radius dielectric particle in the focal plane of 
a highly focused Laguerre- Gaussian beam (reproduced from [33] ). 



frame the stochastic differential equations ([TJ are separated and without loss of generality the treatment can be 
restricted to the x-projection of the system. 

The autocorrelation function (ACF) of the solution to equations (JlJ in each direction reads 

< x(t)x(t + At) >= D^-e- k "\ At ^, (2) 
k x 

where k x is the trap stiffness. Experimentally the trap stiffness is found by fitting the theoretical ACF ^ to the one 
obtained from the measurements. Using the Wiener-Khintchine theorem, the power spectral density (PSD) can now 
be calculated as the Fourier transform of the ACF: 

Pxif) = 2tt2 (p + ff) ' (3) 

where f c — k x /2n r y is the corner frequency. 

A constant and homogeneous external force f ex t,x acting on the probe produces a shift in its equilibrium position 
in the trap. The value of the force can be obtained as: 

fext,x = k x < x(t) >, (4) 

where < x(t) > is the probe mean displacement from the previous equilibrium position. 

The PFM has been applied to measure forces in the range of femto- to pico-Newton in many different fields with 
exciting applications, for example, in biophysics, thermodynamics of small systems, and colloidal physics j4j [6l |7j 
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For a correct use of the PFM, the force field to measure has to be invariable (homogeneous) on the scale of the 
Brownian motion of the trapped probe, i.e. in a range of 10s to 100s of nanometers depending on the trapping 
stiffness. This condition implicates that the force field must be conservative, excluding the possibility of a rotational 
component. However, there are cases where these assumptions are not fulfilled as it is illustrated in Fig. [T] The 
field can vary in the nanometer scale, for example, in the presence of a radiation force field produced by a surface 
plasmon polariton [3D]. It can also be non-conservative in the presence of a rotational force (torque). These effects 
appear, for example, considering the radiation forces exerted on a dielectric particle by a patterned optical near-field 
landscape at an interface decorated with resonant gold nanostructures [3TJ (Fig. [JJa)); the nanoscale trapping that 
can be achieved near a laser-illuminated tip [32] (Fig. |ljb) ) ; the optical forces produced by a beam which carries 
orbital angular momentum [331 (Fig- [ljc)); or in the presence of fluid flows [34"] . 

In this article, we show how to improve the PFM technique in order to be able to deal with these case s. We 
introduce the theory of this enhanced PFM (section [n| . Based on this theoretical understanding, in section III we 



propose a concrete analysis workflow to reconstruct the force field from the experimental time-series of the probe 



position. Finally, in section IV we present experimental results for some particularly important cases, namely the 



case of a conservative or rotational force-field. 



II. THEORY 



In the presence of an external force field f ex t( r (i)), equation ([TJ can be written in the form: 

7*(*) = f(r(t))+V2^7h(t), (5) 
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where the total force acting on the probe f (r(t)) = f e xt(r(t)) — Kr(f) = [f x (r(t)), f y (r(t))] depends on the position 
of the probe itself, but does not vary over time. We reduce our analysis to a bidimensional system, because it is the 
most interesting from the applied point of view. However, our approach can be generalized to the tridimensional case. 
The force f(r(i)) can be expanded in Taylor series up to the first order around an arbitrary point f: 



f(r(i)) 



" /*(?) " 






+ 



d x f x (v) d v f x (r) 
d x f v {v) d v f v {v) 



(r(t)-r) 



(6) 



Jf 



where ff and Jf are the zeroth-order expansion, i.e. the force field value at the point f , and the Jacobian of the force 
field calculated in r, respectively. 

In a PFM the probe particle is optically trapped and, therefore, it diffuses due to Brownian motion in the total 
force field (the sum of the optical one and the one under investigation). If ff ^ 0, the probe experiences a shift in the 
direction of the force. After a time has elapsed, therefore, the particle settles down in a new equilibrium position of 
the total force field, such that ff = 0. As seen in the introduction, the measurement of this shift allows one to evaluate 
the homogeneous force acting on the probe in the standard PFM. Assuming, without loss of generality, r — 0, the 
statistics of the Brownian motion in the surroundings of the equilibrium point can be analyzed in order to reconstruct 
the force field up to its first-order approximation. 



A. Brownian motion near an equilibrium position 

The first order approximation to equation ([5| near a stable force field equilibrium point, f = 0, is: 

r(t) = 1 - 1 3 r(t) + V2Dh(t), (7) 

where r(t) = [x(t), y(t)], h(t) — [h x (t),h y (t)], and J is the Jacobian calculated in the equilibrium point. According to 
the Hclmholtz theorem, under reasonable regularity conditions any force field can be separated into its conservative 
(i.e. irrotational) and non-conservative (i.e. rotational or solcnoidal) components. The two components can be 
immediately identified if the coordinate system is chosen such that d y f x (0) — —d x f y (0). In this case, the Jacobian 
J normalized by the friction coefficient 7 reads: 

7~ 1 Jo = 

where tj> x = k x /-y and <f> y = k y /j, k x = -d x f x (r) and k y = -d y f y (r), and = 7~ 1 <9< / / a; (f ) = -■y~ 1 d lc fy(f). In Q the 
rotational component, which is invariant under a coordinate rotation, is represented by the non-diagonal terms of the 
matrix: f2 is the value of the constant angular velocity of the probe rotation around the z axis due to the presence 
of the rotational force field [33 . The conservative component, instead, is represented by the diagonal terms of the 
Jacobian and is centrally symmetric with respect to the origin. Without loss of generality, we impose that k x > k y , 
i.e. <p x > 4> y . This means that the stiffness of the trapping potential is higher along the x-axis. 
The equilibrium point is stable if: 

rDet(Jo)^0 2 -A(/. 2 + r! 2 >O 

\ Tr (J ) = -2(f) < ' (J) 

where <fi = (4> x + 4> y )/2 and A<p = (<f> x — <j} y )/2. The fundamental condition required to achieve stability is > 0. 
Assuming that this condition is satisfied, the behavior of the optically trapped probe can be explored as a function of 
the parameters f2/</> and A0/</>. The stability diagram is shown in Fig. [2] The standard PFM corresponds to A(j> = 
and = 0. When a rotational term is added, i.e. and A0 = 0, the system remains stable [33] • When there 

are no rotational contributions to the force field (f2 = 0) the equilibrium point becomes unstable as soon as Acf> > 
This implicates that <p y < 0, and therefore the probe is not confined in the y-direction any more. In the presence of a 
rotational component (f2 ^ 0) the stability region becomes larger; the equilibrium point now becomes unstable only 
for A(f> > y/(/) 2 ~n 2 . 

Some examples of possible force fields are presented in Fig. [3j 

When = the probe movement can be separated along two orthogonal directions. As the value of Acf> increases 
the probability density function (PDF) of the probe position becomes more and more elliptical, until for Aip > <fi the 
particle is confined only along the a;-direction and the confinement along the y-direction is lost. 
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FIG. 2: (Color online) Stability diagram. Assuming (f> > 0, the stability of the system is shown as a function of the parameters 
Cl/(j) and A(j>/(j). The regions that satisfy the conditions described in equation Q. The continuos line delimits the stability 
region. The dashed lines represent the A<f> = |0| and A<j> = <f>. The dots represent the parameters that are further investigated 
in Fig. (3), Q and (jj). 
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FIG. 3: (Color online) The arrows show the force field vectors for various values of the parameters A(f>/(f> and Q/cj>. The 
shadowed areas show the probability distribution function (PDF) of the probe position in the corresponding force field. 



If A<f> = 0, the increase in f2 induces a bending of the force field lines and the probe movement along the x- 
and y-directions are not independent any more. For values of fl > cj>, the rotational component of the force field 
becomes dominant over the conservative one. This is particularly clear when Aip ^ 0: the presence of a rotational 
component covers the asymmetry in the conservative one, since the probability density distribution assumes a more 
rotationally-symmctric shape. 



B. Enhanced Photonic Force Microscope 



As we already mentioned in the introduction, the most powerful analysis method is based on the study of the 
correlation functions - or, equivalently, the power spectral density - of the probe position time-series. In this subsection, 
we derive the correlation matrix first in the coordinate system considered in the previous subsection, where the 
conservative and rotational components are separated. We, then, derive the same matrix in a generic coordinate 
system and identify some functions that are independent on the choice of the coordinate system. For completeness, 
we will also present the power spectral density matrix. 
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1. Correlation Matrix 



The correlation matrix of the probe motion near an equilibrium position can be calculated from the solutions of 

r -I T 

(7|, whose eigenvalues are Ai^ = — ± \J A<p 2 — fl 2 and whose eigenvectors are v x 2 = A/p ± \J A<fi 2 — £! 2 
Treating h(t) as a driving function, the solution of ^ is given by: 



where 



W(i) 



r(t) = V2D [ W(t)W- 1 (s)h(s)ds 

J — oo 



n n 

A0+ y/A<j> 2 -fl 2 A(j) - yfA(j) 2 - fl 2 



e Xlt 
e A2t 



(10) 



(11) 



is the Wronskian of the system. 

Since we are assuming r(t) to be a stationary stochastic process, the correlation matrix (r(t + At)r* (t)} can be 
obtained by taking the ensemble average (r(Ai)r* (0)): 



(r(At)r h (0)) = ( 2D / W(At)W~ 1 (t / )h(t , )dt / h h (<")W- lh (t")W h (0)dt A 



where the superscript h indicates the hermitian. Solving this system, we have 



r xx {At) = D- 



-0|At| 



r yy (At) = D- 



,-4>\At\ 



r xy {At) = D 



r yx (At) = D 



e -0|At| n 



e -*|At| n 



fl 2 -a 2 A<J} 2 2 A(j) 

n 2 -Acf) 2 ~ a ~f 

n 2 ~a 2 A<p 2 2 A<p\ 

n 2 -Ac/) 2 +a ~f) 

2 A4> 



C (At)-c 2 ^(l 



2 A0 



C(At) + a z -f 1 



A<t> 

T 

A0\ 



S{\At\) 
S{\At\) 



where 



-S(At) + a 1 -^ (C(Ai) + 5(|At|)) 
-S(Ai) + a 2 ^ (C(At) + 5(|At|)) 



a 



4> 2 + (ft 2 - Acj> 2 ) 



(12) 

(13) 
(14) 
(15) 
(16) 

(17) 



is a dimensionless positive parameter, 

C{t) = 

and 

S(t) = 



cos(\/\A<j) 2 -tt 2 \t) fl 2 >Acj) 2 

i ' n 2 = A0 2 



(18) 



cosh (y|A0 2 -ft 2 |iJ ft 2 < A<t> 2 
sh^^/iA^-n 2 !*) 



sin 

> ^ , 

<pt 

sinh(^/|A0 2 -n 2 |t) 



ft 2 > A^ 2 

n 2 = A0 2 

ft 2 < A</> 2 



(19) 



In Fig. [4] these functions are plotted for different ratios of the force field conservative and rotational components. 
Some cases have already been studied experimentally. For the case A<j> = [33], the ACFs and cross-correlation 
functions (CCFs) are r xx {At) = r yy (At) = De^l A *l cos(nAi)/^ and r xy (At) = -r yx (At) = De^l At l sin(fiAf )/(£, 
respectively. As the rotational component becomes greater than the conservative one (fl > (/>), a first zero appears 
in the ACFs and CCFs and, as ft increases even further, the number of oscillation grows. Eventually, for ft ^> 4> the 
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FIG. 4: (Color online) Auto- and cross-correlation functions for various values of the parameters A<j>/<j> and ip/<f>: r xx (black 
continuous line), r yy (black dotted line), r xy (blue continuous line), and r yx (blue dotted line). 



sinusoidal component becomes dominant. The conservative component manifests itself as an exponential decay of the 
magnitude of the ACFs and CCFs. 

When D, = 0, the movement of the probe along the x- and y-directions becomes independent. The ACFs behave 
as r xx (At) = De-^ At \/(j) x and r yy (At) = £»e~*» |At| /^, while the CCFs are null, r xy (At) = r yx (At) = 0. In Fig. [2] 
this case is represented by the line O = 0. 

When both f2 and A<fi are zero, the ACFs are r xx (At) = r yy (At) = De-^ At \/(/), and the CCFs are null, r xy (At) = 
r yx (At) = 0. The corresponding force field vectors point towards the center and are rotationally symmetric. 

It is also interesting to consider the intermediate cases. In these cases the effective angular frequency that enter the 
expression is given by ^/|A0| 2 — f) 2 . This shows that the difference in the stiffness coefficients along the x- and y-axes 
effectively influences the rotational term, if this is present. A limiting case is when |f2| = A<fi. This case presents a 
kind of resonance between the rotational term and the stiffness difference. However, it is not a dramatic resonance, 
as it is shown by the corresponding force field (Fig. [31. 



2. Correlation matrix in a generic coordinate system 



The expression for the ACFs and CCFs (13 1 to (16) were obtained in a specific coordinate system, where the 
conservative and rotational component of the force field can be easily identified. However, typically the experimentally 
acquired time-series of the probe position required for the calculation of the ACFs and CCFs are given in a different 
coordinate system, rotated with respect to the one considered above. If a rotated coordinate system is introduced, 
such that: 

(20) 







cos 


9 -sin6»" 




X 




Y_ 




sin ( 


? cos 9 




y _ 
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are 


obtained 


as linear combinations of 



r x i x i(At) = (cos#) r xx (At) — cos 9 sin Or xy (At) — svaO cos Or yx ( At) + (sin 9) r yy (At), 
r y / y i(At) = (sin 0) 2 r xx (At) + sin 9 cos 9r xy (At) + cos 9 sin 0r yx (At) + (cos 0) 2 r yy (At), 
(At) = cos sin Or xx ( At) + (cos 9) 2 r xy (At) — (sin 9) 2 r yx (At) — sin cos 9r yy ( At), 



x'y 



y'x 



(At) = sin 9 cos 9r xx (At) - (sin 9) 2 r xy ( At) + (cos t 



' yx 



(At) 



cos sin t/r 



9r yy (At), 



(21) 
(22) 
(23) 
(24) 



which in general depend on 
ACFs are invariant: 



x'y 



(At) 



However, it is remarkable that the difference of the two CCFs and the sum of the 

g-0|At| Q 

S(At), 

n 2 -a 2 A<j) 2 \ „ /A ^ , 2 Aah 2 



' y'x 



(At) = 2D- 



' x'x' 

(At) +r y , y ,(At) = 2D 



,-4>\M\ 



n 2 ~ Ac/) 2 



C(At) 



-S(\At\) 



(25) 
(26) 
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FIG. 6: (Color online) Functions which depend on the orientation of the coordinate system for various values of the parameters 
as a function of the angle with respect to the coordinate system chosen in the section]!!] (a) r x 'x' — r y'y' an <A (b) r x i y i + r y > x >. 
The thicker lines show the values for At = 0. 



These functions are presented in Fig. [5] These functions are very similar to the ones presented in Fig. [4] however, 
the latter depend on the coordinate system choice. 

Other two combinations of (21 l-p4]), which are also useful for the analysis of the experimental data, namely the 
sum of the CCFs and the difference of the ACFs, depend on the choice of the reference frame: 

e -<f>\At\ AJ, / Q \ 

r x , y ,{M) + r y>x> {M) = 2D^^a 2 ^ (C(At) + S(\At\)) (-cos(28) - sm(26)j , (27) 

-0|At| A A / Q \ 

r x , x ,{M)-r y , y ,{&t) = -2D — - — a 2 -f- (C(At) + S(\At\)) - sin (29) + cos (20) . (28) 



Their plots are shown in Fig. [6] In particular, when they are evaluated for A< = 0, they deliver information on the 
orientation of the coordinate system. 
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3. Power Spectral Density Matrix 

In the frequency domain the equation ^ is given by: 

;2tt/R(/) = J R(/) + V2DK(f), 



(29) 



and its solution is R(/) = \/2D (i2nfl2 — Jo) 1 H(/), where I2 is the 2D unit matrix, and the corresponding PSD 
matrix: 



P(/) = R R h = 



2D 



b x + i2nf) {4> v + i2nf) + ijj 



•>\> 



^+47r 2 / 2 + ^ 2 i>[cf> x -<t> y -i4nf] 
i> [4> x -<j>y + i4nf\ eg + A^f 2 + ^ 



(30) 



where the property H(/) • H h (/) = I 2 has been used. We notice that the PSD matrix could have been obtained as 
Fourier-transfor of the correlation matrix ( Wiencr-Khintchine theorem) . 



III. EXPERIMENTAL CONSIDERATIONS 



In this section we propose a concrete analysis workflow to reconstruct the force field from the experimental time- 
series of the probe position. 

Experimentally the probe position time-series is the only available information to reconstruct the force field. Typ- 
ically these data are obtained in an arbitrary coordinate system. These time-series need to be statistically analyzed 
in order to reconstruct all the parameters of the force field, i.e. cf>, A(f>, and fl, and the orientation of the coordinate 
system. The detailed procedure to retrieve all this information from the experimental data is presented in this section. 

Let us suppose to have the probe position time-series in a generic coordinate system r'(t) = [x'(t),y'(t)], First, 
we evaluate the parameters <f>, A(f>, and fl. Then, we transform the coordinate system to the one presented in the 
section [IT] where the conservative and rotational components are separated. Finally, we reconstruct the total force 
field. Eventually, the trapping force field may be subtracted to retrieve the external force field under investigation. 

In order to illustrate this method we proceed to analyze some numerically simulated data. The main steps of this 
analysis are presented in Fig. [7] In Fig. [7](a) the PDF is shown for the case of a probe in a force field with the 
following parameters: (f> = 37.1 s _1 , A(f> = 9.3 s -1 (corresponding to k x — A3.25pN//im and k y — 26.25 pN/fim), 
= 0, and 9 — 30°. The PDF is ellipsoidal due to the difference of the stiffness along two orthogonal directions. In 
Fig. Pnb) the PDF for a force field with the same cf> and Atfi but with fl = 37.1 s _1 is presented. The presence of the 
rotational component in the force field produces two main effects. First, the PDF is more rotationally-symmetric and 
its main axes undergo a further rotation. Secondly, as we show below, the CCF is not null (Fig. Rd)). 



A. Evaluation of the parameters (f>, A(f>, and fl 



In order to evaluate the parameters of the force field, 
coordinate system where the experiments have been done, 



A<p, and fl, we calculate the correlation matrix in the 



r'(At)r' h (0)W 



.(At) r xV (Ai) 



x'(At) 



y'y' 



(At) 



(31) 



Then we calculate the CCF difference (25 1, r x > y i(At) 



y'x 



(At). As we showed in section II, this function is invariant 
with respect to the choice of the reference system, and it is different form zero only if O ^ 0. The results are shown 
in Fig. [7jc) and [7jd) for the cases of the data shown in Fig. [7ja) and [7jb) respectively. The three aforementioned 
parameters can be found by fitting the experimental data to the theoretical shape of this function. In particular, the 
exponential decay of the function is related to the <j> parameter; the period of the superimposed oscillations is related 
to the effective angular frequency yJ\A(f) 2 — f2 2 |; and the sign of the slope in At = gives the sign of fl. 

When O = 0, the CCF difference p5| is null (Fig. (7^c)), it can not be used to find the two remaining parameters. 
For fl = 0, the other invariant function, the ACF sum (equation ( 26 1 ) , is given by 



'x'x' 

(At)+r y , y ,(At) = 2D- 



-0|At| 



A<t> 2 



C(At) 



(32) 



and A<j) can be evaluated by fitting the data to (32). The function (26 1 can be used for the fitting of the three 



parameters but can not give information on the sign of fl, which must be retrieved from the sign of the slope at 
At = of the CCF difference. 
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FIG. 7: (Color online) (a-b) Probability density function for a Brownian particle under the influence of the force-field shown 
in the inset; in (a) the force-field is purely conservative, while in (b) it has a rotational component, (c-d) Invariant function, 
f'x'x 1 + r y'y' (black line) and r x i y i — r v i x i (blue line), (e-f) Force fields obtained form the numerically simulated data. 



B. Coordinate system transformation 



Although the values of the parameters <p, A</>, and f2 are now known, the directions of the force vectors are still 
missing. In order to retrieve the orientation of the experimental coordinate system, we now use the orientation 
dependent functions (27 1 and (28). The best choice is to evaluate the two functions for At = 0, because the signal- 



to-noise ratio is highest at this point: 



r* v (0) + Vx'(0) = 2 D i^f (? cos (29) - sin (28)) 
TV^O) - r yV (0) = -2Df^ (f sin (20) + cos (26) 



(33) 



The solution of this system delivers the value of the rotation angle 9. If A(f> — 0, (33 1 is undetermined as a consequence 
of the PDF radial symmetry. In this case any orientation can be used. If f2 = 0, the orientation of the coordinate 
system coincide with the axis of the PDF ellipsoid and, although (33 1 can still be used, the Principal Component 



Analysis (PC A) algorithm is a convenient means to determine their directions. 



C. Reconstruction of the force field 



Now everything is ready to reconstruct the unknown force field acting on the probe around the equilibrium position 
in an area comparable with the mean square displacement of the probe. From the values of cj> and A<fi, the conservative 
forces acting on the probe result f c (x, y) = — (k x xe x + k y ye y ) and, from the values of f2, the rotational force is 
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f r (x, y) — ft (ye x — xe y ). The total force field is 



f (or, y) = f c (x, y) + f r (x, y) = (-k x x + fly) e x - (k y y + fix) e y x 



(34) 



in the rotated coordinate system. Now the rotation (20 1 can be used in order to have the force field in the experimental 
coordinate system. The unknown component can be easily reconstructed by subtraction of the know ones, such as 
the optical field. 



IV. EXPERIMENTAL RESULTS 



For an experimental verification of our conclusions, we analyze the Brownian motion of an optically trapped 
polystyrene sphere in the presence of an external force field generated by a fluid flow j34j . A schematic of the setup 
is presented in Fig. [8] 

An optical trap is generated by a CW 633 nm beam at the focal plane of a lOOx 1.3AM objective lens inside 
a chamber. The chamber is prepared using two cover slips separated by a 50 fim spacer and filled with a solution 
containing polystyrene spheres (radius R = 0.5 /jm). The forward scattered light from the trapped sphere is collimated 
by a 50 x objective onto a quadrant photodiode (QPD). The trap force constant can be adjusted by changing the 
intensity of the laser beam. 

The fluid flow that produces the external force field was generated using solid spheres made of a birefringent material 
(Calcium Vaterite Crystals (CVC) spheres, radius R = 1.5±0.2/um [35]), which can be made spin due to the transfer 
of orbital angular momentum of light. They are ail-optically controlled, i.e. their position can be controlled by an 
optical trap and their spinning state can be controlled through the polarization state of the light. In our experimental 
realization up to four CVC spheres were optically trapped in water and put into rotation using four steerable 1064 nm 
beams from a Nd:YAG laser with controllable polarization - to control the direction of the rotation - and power - to 
control the rotation rate. 



Conservative force field 



In order to produce a conservative force field, two CVC were placed as shown in Fig. (9ja), which should theoretically 
produce the force field presented in Fig. |9jc). 

In Fig. 10 'a), the invariant functions, r x / x i + r y i y i (black line) and r x i y > — r y i x > (blue line), and respective fitting 



to the theoretical shapes are presented. The CCF difference tells us that ft = in this case, while the fitting to the 
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FIG. 9: (Color online) Experimental configuration with two (a) and four (b) spinning beads and respective force fields (c) and 
(d). 
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FIG. 10: (Color online) (a) Invariant functions, 7V X ' + r y 
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(black line) and r x 'y< 



(blue line), and respective fitting to 



the theoretical shape (dotted lines), (b) Experimental probability density function and estimated force field. 

ACF sum tells us the values of = 18 s^ 1 and A(j> = 6 s _1 . The value of the rotation of the coordinate system in this 
case is 32°. 

The total force field can now be recontructed: k x — 225 fN//j,m and k y = 112 fN/ \im. This force field is presented in 
Fig. 10 'b). We can now retrieve the hydrodynamic force field by subtracting the optical force field (k opt — 185 fN/fim 
approximatively constant in all directions), that can be measured in absence of rotation of the spinning beads. 



B. Rotational force field 

In order to produce a rotational force field, four CVC were placed as shown in Fig. |9jb), which should theoretically 
produce the force field presented in Fig. |9jd). 

In Fig. 11 a), the invariant functions, r x i x > + r y i y i (black line) and r x > y i — r y i x i (blue line), and respective fitting to 
the theoretical shapes are presented. Now the CCF difference is not null any more and tehrefore it can be used to fit 
the three parameters: 0=11 s _1 , A<fi w 0, and f2 = 5 rads" 1 . We can notice that the ACF sum can be used for this 
purpose as well; however, we have to remark that using the latter the sign of stays undetermined. The small value 
of A</> implicates that the rotation of the coordinate system is not crucial. 



The total force field can now be recontructed: k x ss k y — 100 fN/ fim. This force field is presented in Fig. 11 'b). We 
can now retrieve the hydrodynamic force field by subtracting the optical force field (k opt — 78 fN/fim approximatively 
constant in all directions), that can be measured in absence of rotation of the spinning beads. 
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(black line) and 



(blue line), and respective fitting to 



the theoretical shape (dotted lines), (b) ) Experimental probability density function and estimated fore field. 



V. CONCLUSION 



We have shown how the PFM can be applied to the detection of locally non-homogeneous force fields. This has 
been achieved by analyzing the ACFs and CCFs of the probe position time-series. We believe that this technique can 
gain new insights into micro- and molecular-scale phenomena. In these cases the presence of the Brownian motion is 
intrinsic and has can not be disregarded. Therefore this technique permits one to take advantage to the Brownian 
fluctuations of the probe in order to explore the force field present in its surroundings. 

One of the most remarkable advantages of the technique we propose is that it can be implemented in all existing 
PFM-setups and even on data acquired in the past. Indeed, it does not require changes to be made in the physical 
setup, but only to analyze the data in a deeper way. 
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